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Castaing's instability in a trapped ultra-cold gas 
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Abstract. We consider a trapped ultra-cold gas of (non-condensed) bosons with two internal states (de- 
scribed by a pseudo spin) and study the stability of a longitudinal pseudo spin polarization gradient. For 
this purpose, we numerically solve a kinetic equation corresponding to a situation close to the experiment 
at JILA jjj]. It shows the presence of Castaing's instability of transverse spin polarization fluctuations at 
long wavelengths. This phenomenon could be used to create spontaneous transverse spin waves. 

PACS. 03.75.Fi Phase coherent atomic ensembles; quantum condensation phenomena - 51.10.+y Kinetic 
and transport theory of gases - 75.30.Ds Spin waves 



1 Introduction 

Recent experiments |l| with a trapped ultra-cold (non 
, condensed) Bose gas with two internal states have shown 
the existence of interesting relative population dynamics. 
| Theoretical studies [|| || making use of a pseudo spin de- 
: scription for the internal states have explained this phe- 
nomenon as pseudo spin oscillations due to the "identi- 
cal spin rotation effect" (ISRE) which appears when 
the temperature is low enough for the binary collisions 
to be in the quantum regime. Alternatively, this mecha- 
nism can be understood as a "spin mean- field" ||. In a 
polarized system, this implies the existence of low energy 
excitations of the transverse spin polarization, named spin 
■ waves. The prediction of spin waves in dilute gases §§ 
as well as their observations in H| and helium goes 
. back to the 1980's. Shortly later, Castaing || showed that 
a strong gradient of longitudinal spin polarization is un- 
stable with respect to transverse fluctuations. His study 
focused on homogeneous polarized 3 He gas and assumed 
that the spin waves were in the hydrodynamic regime. The 
purpose of this article is to provide a quantitative study 
of the existence of Castaing's instability in an inhomo- 
geneous system relevant to the experimental situation at 
JILA [Q. In this experiment, neither the hydrodynamic 
nor the collisionless regime for the spin oscillations are 
reached. This work is motivated by the recent contribu- 
tion of Kuklov and Meyerovich , who were the first to 
suggest and study the existence of Castaing's instability 
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in this context. For completeness, one should be aware of 
the numerical studies of Castaing's instability in polar- 
ized Fermi fluids, which are based on the solution of the 
Leggett equations rather than the full kinetic equation 



2 Kinetic equation; Leggett equations 

The physical situation we consider is close to that of the 
experiment at JILA Q: 87 Rb atoms (bosons) with two hy- 
perfine states of interest (denoted by 1 and 2) are confined 
in an axially symmetric magnetic trap elongated in the Ox 
direction. The temperature T is about twice the critical 
temperature for Bose-Einstein condensation, so that the 
gas is non-degenerate (Boltzmann gas). However, the de 
Broglic thermal wavelength is much larger than the scat- 
tering length so that the collisions occur in the quantum 
regime. It is convenient to consider the pseudo spin asso- 
ciated with the two hyperfine states 1 and 2 (the basis in 
spin space is denoted by {ej_,i; ej.,2! e ||})- Initially the spin 
polarization is longitudinal (i.e. along the "effective exter- 
nal magnetic field" ) and has a strong spatial gradient. For 
example, on the left (resp. right) of the trap center, the 
cloud of atoms is mostly in state 1 (resp. 2). This might 
be achieved, for example, by preparing a cloud of atoms 
in state 1 and a cloud of atoms in state 2 separated by a 
sharp optical potential barrier at the center of the mag- 
netic trap. After removal of the optical barrier, as the two 
clouds mix, a strong longitudinal spin polarization gradi- 
ent appears in the region of overlap. If it is strong enough, 
one should be able to observe the appearance of a large 
transverse component of the spin polarization as a result 
of Castaing's instability. 

To study this system^we write an effective one dimen- 
sional kinetic equation pi in terms of a local density in 
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phase space f(x,p, t) and (pseudo) spin density M.(x,p, t), 
where x is the position, p the momentum, and t the time. 
The one dimensional description for an elongated system 
is justified by the time scale separation associated with ra- 
dial and axial motions characterized by frequencies Lo ra d 
and Uax respectively. The time scales differ by more than 
an order of magnitude leading to an effective dynamical 
averaging over radial dynamics as described in Ref. (U || . 
We define the local density n and spin polarization m by 



(1) 



n(x,t) = J dpf(x,p,t) 
m (s, t ) = /*«(*,,,«) 



The local density of atoms in state 1,2 is then given by 
n i.2 = (n^ mii)/2, defining m\\ as the longitudinal com- 
ponent of the spin polarization. We also denote g+j = 
4irh 2 aij /m as the coupling constants associated to the dif- 
ferent scattering lengths a«j where i and j — 1,2 (a 2 i = 
ai2 and to is the mass of the particles). 

We now review the different terms entering the kinetic 
equation (for details, see Ref. 0). The force acting sim- 
ilarly on both internal states contains three terms: the 
magnetic trap, the mean-field, and the Stern-Gerlach force 
(associated to a gradient of "effective external magnetic 
field" f2, see below). The magnetic trap force is dominant, 
so the associated potential energy is given by 



Vl(x) + V 2 (x) 1 2 2 

= -mw„j 

2 2 ax 



(2) 



where V± and V 2 are the magnetic trapping potentials act- 
ing on states 1 and 2. The differential force can be de- 
scribed in the pseudo spin picture by an "effective exter- 
nal magnetic field" fl which contains two contributions: 
one due to the differential Zeeman effect and one due to 
the differential mean field, so that 

M2{x) = [V 2 {x) - V x (x) + ( 522 - ffn)n(i)/2]e|| (3) 

where en is a unit vector in the longitudinal direction. Fol- 
lowing Ref. 0], we have assumed that 2gi2 ~ 511 + 1722 for 
simplicity. The average value over the sample of the "ef- 
fective external magnetic field" is removed by going to a 
uniformly rotating frame (Larmor frame) . It is also crucial 
to consider the ISRE which manifests itself as a "molec- 
ular field" or "spin mean field" 512 m(x,t)/2, which adds 
to the "effective external magnetic field" (m. Finally, the 
collision integral is treated in a relaxation-time approx- 
imation with a time r of the order of the average time 
between collisions. The kinetic equation is written as 

d t f+^d x f-mu> 2 ax xd p f~ -{f-D/r (4) 

d t M + — d x M - muj 2 x d p M -(12+ ^^) x M 

~ -(M-M eq )/ T (5) 

where f eq and M ei? are local equilibrium phase space den- 
sities. 



In the hydrodynamic regime close to local equilibrium, 
the above equations reduce to the counterpart of the Leggett 
equations Jll| in the case of a Boltzmann gas || , namely 

dfjn + d x j = fl x m 



a . (rt , 512m . k B T 
d t j -{S2+ -~r-) x J + — d x m ■ 



m 
j 

r 



(0) 



where t) is the spin polarization current along Ox. 
These equations are also valid in the coUisionless regime 
close to global equilibrium Q. We define /x = gi 2 n(0)T/2H 
as the dimensionless parameter which characterizes the 
strength of the ISRE and denote D = ksTT/m as the 
spin diffusion coefficient. 



3 Castaing's instability 

Studying polarized 3 He gas, Castaing || noticed that a 
sufficiently strong longitudinal spin polarization gradient 
is unstable against transverse long wavelength fluctua- 
tions. This result was obtained by analyzing a uniform 
gas (uJax — 0) with a uniform precession frequency (so 
that fl = in the Larmor frame) and by assuming a time 
independent spin gradient. Using the Leggett equations 
(||), Castaing studied the stability of this system with re- 
spect to the transverse fluctuations. He assumed a small 
plane wave perturbation around a slightly non-uniform 
stationary solution such that 

m(x,t) = m|}(aj)e|| + Sme^ kx -^ 
j(x, t) = -Dd x m\{x) H + fle****-"*) (7) 

For circularly polarized transverse spin waves in the pres- 
ence of a strong gradient of longitudinal spin polarization, 
one obtains the following spectrum for the hydrodynamic 
regime (d t j a 0) 



D 



1 + (fim?Jn) 2 



^/imjj/n - ij {k 2 - fik(d x m°)/i 



(8) 

Here, a mode with wave vector k is unstable if the mode 
frequency lo has a positive imaginary part u>i, which leads 
to instability when 



k 2 < nk(d x m?A/n 



(9) 



The goal of this paper is to study Castaing's instability 
in a trapped ultra-cold gas. Using a numerical simulation, 
described in the next section, we show that an instability 
does occur. The situation we study is richer than the one 
considered in Castaing's original work in several aspects: 

(i) with the full kinetic equation (instead of Leggett equa- 
tions) we can describe spin oscillations of large amplitude, 

(ii) we are not limited to the hydrodynamic regime (for 
example, regimes between hydrodynamic and coUisionless 
are included), (iii) the longitudinal spin polarization gra- 
dient is not assumed to be constant, and (iv) we can also 
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Fig. 1. Initial density of atoms in state 1 (solid line) and 
2 (dashed line) giving a strong gradient of longitudinal spin 
polarization near the center of the trap (d x m®/n(0) — 2/xt)- 
The total density is also plotted (dotted line). 

include a possible gradient of the precession frequency (i.e. 
the "effective external magnetic field" is position depen- 
dent). 



4 Numerical simulation 

The kinetic equation (^,||) is solved numerically by prop- 
agating in time the initial distribution in a discretized 
phase-space using the Lax-Wendroff method (see e.g. fli^l). 
We assume that the density in phase space is initially 
at equilibrium / = f eq cx exp(— x 2 /2x T — p 2 /2p T ), with 
xt — /mu! 2 x and pr = y/mkBT. As the equilibrium 

distribution does not evolve in time under the action of 
(|J) , the local density is constant and we concentrate on the 
evolution of M described by the kinetic equation (|J) . The 
initial spin density distribution M is taken as the product 
of the Maxwell-Boltzmann equilibrium distribution along 
the longitudinal axis and a function which equals —1 to 
the left of the trap center (x < —xt), +1 to the right 
(x > xt) and has a constant positive slope through the 
center (— xt < x < xt)- The resulting initial density of 
atoms in state 1 and 2 is shown in Figure |l|. We introduce a 
small initial transverse perturbation in order to start the 
instability. Without loss of generality, the perturbation 
can be written as a plane wave (with spatial frequency of 
order 1/xt) multiplied by the Maxwell-Boltzmann equi- 
librium distribution and divided by a number M S> 1 
(typically between 10 3 and 10 6 ): 

M ±A {x,p) = H(x lP )cos(Trx/x T )/JV (10) 

Parameters used in the simulation are taken from Ref. 

The axial trapping frequency is uo ax /2TT = 7 Hz and 
the time between collisions r ~ 10 ms, so that Lo ax T ~ 
0.5. The "effective external magnetic field" is taken to 
be an inverted Gaussian of depth 5f2 (\8fi\/u> ax is varied 



Fig. 2. Time evolution of the (normalized) transverse spin 
polarization mi,i(x,t)/mi,i(a;,0) at different positions in the 
trap: x/xt = (solid line); 5/25 (dashed line); 10/25 (dotted 
line); 15/25 (dashed-dotted line). For this simulation /i ~ 22, 
uj ax T ~ 0.6, 8Q/uj ax = 2, d x m\/n(Q) = 2/x T and TV = 10 3 . 



between and 2) and half- width xt- The density at the 
center of the trap is n(0) = 1.8 x 10 19 m -3 . The initial spin 
polarization gradient near the center of the trap is varied 
and typically |9 x mii|/n(0) ~ 2/xt- For 87 Rb atoms, with 
a\2 = 5.2 x 10~ 9 m and T = 0.6 fiK, the ISRE parameter is 
/i ~ 5 . With these parameters, an instability can already 
be observed in our simulations but is probably too weak 
to be detected experimentally. To enhance the effect, we 
take a scattering length 5 times smaller (ai 2 — * ai 2 /5) and 
an axial trap frequency 20 times smaller (uj ax — > tv ax /20), 
keeping a constant density at the center of the trap. We 
discuss these matters further in the next section. 

As the phase space is modeled by a discrete grid with 
finite size, it is important to distinguish between a physical 
instability and a numerical one. This can be easily done 
as only the latter will depend on the grid spacing. One 
can always get rid of a numerical instability by choosing 
a sufficiently tight grid. 

Results of the simulation are shown in Figures |^ and |^. 
The time evolution of the transverse spin polarization is 
plotted in Figure || for different positions near the center 
of the trap, where the longitudinal spin polarization gra- 
dient is most pronounced. The instability is clearly visible, 
owing to a large enhancement of the transverse spin polar- 
ization by a large factor comparable to Af. Figure || shows 
the logarithm of the time evolution of the (normalized) 
transverse spin polarization at the center of the trap. In 
this representation, one can clearly distinguish the expo- 
nential rise of the envelope (after a short delay of order of 
the time between collisions) and its subsequent exponen- 
tial decay. 
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Fig. 3. Time evolution of the logarithm of the (normalized) 
transverse spin polarization mx,i(0, t)/m± t i(0, 0) at the center 
of the trap. The two lines are plotted to visualize the exponen- 
tial rise of the instability (full line) and then the exponential 
decay of the spin wave (dashed line) . The parameters have the 
same values than for the simulation of Fig. |^. 



5 Discussion 

The initial rise of the transverse spin polarization shows 
that an instability indeed occurs. As already stated, the re- 
sults presented here were obtained for a scattering length 
a\2 and trap frequency uj ax that were smaller than in the 
current experiment at JILA JjJ (by factors of 5 and 20 
respectively). This was done to enhance the effect of the 
"spin mean-field" (to favor the instability) by increasing 
the ISRE parameter (JR. As will be seen below, we have 
to keep WaxT ^ 1, which implies scaling the axial trap 
frequency with the square of the scattering length since 
t -1 oc a\ 2 . This procedure for enhancing the ISRE in- 
sures that the gas remains non-degenerate. 

The phenomenon observed in Fig. || and ^ may be bro- 
ken down into the following four steps: 

I. ) During a time of the order of the time between 
collisions r ~ 250 ms a hydrodynamic description is not 
valid (the spin polarization current, which is a fast vari- 
able, has not reached its stationary value yet). Figures || 
and show that the transverse spin polarization does not 
evolve significantly. 

II. ) Then, as shown by the solid line in Figure ^, the 
envelope of the transverse spin polarization rises expo- 
nentially (the coefficient in the exponential luj is almost 
constant in time). This is a characteristic of an insta- 
bility. The value of u>i (as compared to the formula ob- 
tained by Castaing, see equation (||)) is modified by the 
presence of the trap . An estimate of the time needed 
for the instability to develop (see Ref. |J) is given by 
t 



in si 



L^/D ~ \jur ax T ~ 1 s (where L m is the char- 



1 This could be done experimentally by using a Feshbach 
resonance to decrease the scattering length, as /i ~ A-r/ai2 
where At is the de Broglie thermal wavelength. 



acteristic size of the longitudinal spin polarization gradi- 
ent; L m ~ xt in a typical simulation). It reproduces cor- 
rectly the order of magnitude of the observed maximum 
of transverse spin polarization in the simulation. The in- 
stability can develop only if the spin polarization current 
has reached its stationary value, so that t < ti nst which 
implies u> ax T < 1. 

III. ) The imaginary part of the frequency u>i varies 
slowly in time because the gradient of longitudinal spin 
polarization is not constant but decays as a result of both 
spin diffusion and the presence of the trap. The change of 
sign of u>i marks the end of the exponential grow. Once 
loi is negative, the transverse spin polarization decays as 
an ordinary damped spin wave. 

IV. ) The longitudinal spin polarization gradient finally 
decays on a time scale of the order of the diffusion time 
tdiff ~ L 2 /D ~ 4 s where L, is the size of the cloud 
(L/xt ~ 2). As emphasized by Kuklov and Meyerovich 
P, if L > L m , the instability develops faster than the 
longitudinal spin polarization relaxes, in accordance with 
the results of our simulation. Once the longitudinal spin 
polarization is zero, loi ceases to evolve in time, and conse- 
quently the envelope of the transverse spin wave decreases 
exponentially (dashed line in Figure |^) . This happens for 
t > tdiff - The frequency of the spin wave is not completely 
determined by the real part of the mode frequency in equa- 
tion (||) , as the gradient of external precession frequency, 
dfi, and the presence of the trap are not taken into ac- 
count. Actually, its order of magnitude (at the center of 
the trap) is given by the external precession frequency 
5(1 = 2cu ax ~ 2ir x 0.7 Hz. We have checked that when 
§Q/u) ax = 0, the transverse spin polarization at the center 
of the trap does not oscillate. 

The criterion for an instability (see equation (|j|)) is 
qualitatively verified at the center of the trap if one takes 
into account the fact that wave vectors k are limited by 
the presence of the trap. This forces k > 2tt/L and the 
criterion becomes 



\d x m°\ 



> 



2n 



TV 

XT 



(11) 



Using an initial longitudinal spin polarization gradient 
d x m®/n(0) — 2/ xt and L ~ 2xt implies that the ISRE 
parameter [i should be larger than ~ 2. In the simulation, 
we found an instability threshold at fi ~ 4. Obtaining a 
significant instability requires a much larger value of the 
ISRE parameter. 

We now discuss the relevance of this criterion for one of 
the experiments done at JILA (sec the first article of Ref. 
0), where spatial separation of the two internal states 
was observed as a result of an initial 7r/2 rf pulse. The 
maximum longitudinal spin polarization occured at the 
maximum of spin state separation and can be estimated 
as |<9 x mi| |/n(0) ~ 1/xt (at x ~ ±xt/2). Since the value 
of the ISRE parameter is fj, ~ 5, our numerically obtained 
criterion (fi\d x m^\/n(0) > 8/xt) predicts no instability. 
Whether it is experimentally feasible in practice to reach 
transient total (or nearly total) separation of the two spin 
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states, so that, when re-mixing, the longitudinal spin po- 
larization gradient would be strong enough for Castaing's 
instability to develop, is not clear. A very strong initial 
gradient of longitudinal spin polarization with fj, ~ 5 may 
not be enough to start the instability because the spin 
diffusion is very efficient in decreasing a strong gradient. 
One would have to maintain a strong longitudinal spin 
polarization gradient in order to create an instability with 
the current experimental value of ji. 



6 Conclusion 

Our numerical simulation confirms the possibility of ob- 
serving Castaing's instability in a trapped ultra-cold gas 
with two internal states, as proposed by Kuklov and Me- 
yerovich j^] . We solved a one dimensional kinetic equation 
numerically, and were therefore able to include effects that 
are beyond the usual treatment of Castaing's instability in 
terms of a small amplitude hydrodynamic description. We 
argue that Castaing's instability was probably not relevant 
for the previous experiments done at JILA jjj as both the 
ISRE parameter /i and the longitudinal spin polarization 
gradient were too small. This does not preclude obser- 
vation of the instability in future experiments if relevant 
parameters like the trapping frequency and the scatter- 
ing length are chosen appropriately. We suggest the use 
of Castaing's instability as a way of creating spontaneous 
transverse spin waves as a result of a strong initial longitu- 
dinal spin polarization gradient. Our calculations are also 
valid for a non-degenerate gas of fermions (see Ref. ||), 
where similar effects could be observed, in a case where 
gn — 922 — and the "spin mean-field" changes sign. 
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